Baryon number conservation and the cumulants of the net proton distribution 
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We discuss the modification of tlie cumulants of the net baryon and net proton distributions 
due to the global conservation of baryon number in heavy-ion collisions. Corresponding probability 
distributions and their cumulants are derived analytically. We show that the conservation of baryon 
number results in a substantial decrease of higher order cumulants. Based on our studies, we propose 
an observable that is insensitive to the modifications due to baryon number conservation. 
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I. INTRODUCTION 

The phase structure of the strong interactions, Quan- 
tum Chromodynamics (QCD), has been studied theoret- 
ically and experimentally for many years. 

On the theory side, the problem has been investigated 
within various models and addressed systematically in 
first-principles lattice QCD (LQCD) calculations. Due 
to the sign problem, rigorous results from LQCD on the 
phase structure are presently available only at vanishing 
net baryon density, where it has been established that 
QCD exhibits an analytic crossover transition Jl [with a 
pseudo-critical temperature of Tc ~ 160 MeV [2, 3- 

Exploratory LQCD studies at finite net baryon den- 
sity 043' ^s well as model results (see, e.g. Ref. Q), 
provide some indications for the existence of a critical 
end point (CEP) of a first-order phase transition at fi- 
nite temperature and density. 

In the laboratory, the properties of hot and dense QCD 
matter have been studied in heavy-ion collisions. Since 
the location of the CEP and the first-order coexistence 
region is not really known from theory, a search of this 
phase structure requires the exploration of the entire ex- 
perimentally accessible phase diagram. This is achieved 
by varying the beam energy and monitoring observables 
for possible non-monotonic behavior. A first set of exper- 
iments of this nature has been carried out by the NA49 
collaboration 9, 10] at the CERN SPS and is now pur- 
sued at somewhat higher energy in the recent beam en- 
ergy scan program at RHIC ^llj . Both the CEP and the 
first-order phase transition are associated with character- 
istic fluctuations - long range ones for the second-order 
transition at the CEP, and possible s pin odal instabilities 
in case of a first-order transition [12|, [ij] . The first mea- 
surements by the NA49 collaboration have concentrated 
on the variances, or second-order cumulants, of various 
particle ratios and the transverse momentum as was pro- 
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posed in Ref. [l^. Their measurements showed hardly 
any deviation from the expected Poisson fluctuations of 
a Hadron Resonance Gas (HRG). 

Recently it has been realized that higher order cumu- 
lants, especially of the baryon number, which serves as 
an order parameter of the QCD phase transition at finite 
density, would be more sensitive to the fluctuations as- 
sociated with the second order transition, including the 
CEP [l^- First, because higher order cumulants scale 
with higher powers of the correlation length, a finite and 
limited increase of the correlation length as a result of 
critical slowing down may still be visible. Second, since 
the baryon number is a conserved quantity its fluctu- 
ations are less modified by the final state interaction in 
the hadronic phase 'lq| . Furthermore, model calculations 
found that the kurtosis, i.e., the ratio of the fourth- over 
the second-order cumulant, is negative at high baryon 
chemical potential, close to the chiral crossover line above 
the CEP 17]. In Ref. [l^, the lines, where the kurtosis 
changes sign, were obtained from universality arguments 
in the critical region of the CEP. These theoretical pre- 
dictions suggest a non-monotonic behavior of the kurtosis 
of the baryon number or electric charge distributions as 
a function of the collision energy if chemical freeze-out 
happens close to the CEP. 

In addition, higher order cumulants also provide sensi- 
tive information about the crossover at vanishing baryon 
density. As shown in Ref. [l9|, the sixth-order cumulant 
of the net baryon distribution is negative close to the 
crossover line. Preliminary LQCD results 20] also indi- 
cate negative values of the sixth-order cumulant close to 
the crossover temperature at zero baryon chemical po- 
tential. Thus, by measuring the sixth-order cumulants, 
one may relate the crossover and the freeze-out line at 
non zero chemical potential and provide information on 
an approximate position of the possible CEP. 

When comparing theoretical (model or LQCD) pre- 
dictions to the experimental data, one has to keep in 
mind that the singular behavior of fluctuations is pre- 
dicted in the grand canonical formulation of thermody- 
namics where conservation laws are imposed only on the 
average. Consequently, to address the same physics ex- 
perimentally, one is required to approximately achieve 



conditions of the grand canonical ensemble, i.e., to study 
fluctuations in a restricted phase space ^. This might 
be done by performing appropriate cuts in the rapid- 
ity and/or transverse momentum of detected particles. 
Clearly, the smaller the fraction of observed particles the 
smaller is the effect of global baryon/charge conservation. 
This was also demonstrated in Ref. ^V\ using the UrQMD 
model. However, care has to be taken that these cuts do 
not destroy the underlying correlations responsible for 
the physics one tries to access. Therefore, a fine balance 
between the need to suppress the effects of conservation 
laws and the requirement to preserve the dynamical cor- 



relations has to be found [16|. The subsequent studies. 



we believe, will help in achieving this nontrivial task. 

In this paper, we will explore to what extent global 
baryon conservation modifies cumulants of net baryon 
and net proton distributions as well as their ratios. To 
this end we consider a system where the only correlations 
are due to global baryon number conservation. There- 
fore, we start with Poisson statistics for baryons and 
anti-baryons and, subsequently, enforce baryon number 
conservation. Since the Hadron Resonance Gas (HRG) 
model in the classical or Boltzmann approximation is also 
governed by Poisson statistics for the baryon and anti- 
baryons, our results may be directly applied to the HRG, 
which is commonly used as a theoretical baseline for the 
analysis of heavy-ion collisions ^. In other words, our ap- 
proach is equivalent to a treatment of the HRG in the so 
called canonical ensemble [22] with respect to the baryon 
number. A canonical treatment of the HRG with respect 
to strangeness has been reported in the literature, e.g., in 
Refs. [22, [231 . The effect of a globally conserved charge on 
the variance of the charge distribution has been studied in 
the same framework in Refs. |24|,[25|. Here we will extend 
these studies to higher order cumulants, where the effects 
due to global baryon number conservation are expected 
to be stronger. Therefore, our results may be considered 
as an improved Hadron Resonance Gas prediction for the 
baryon number cumulants, and thus provide a baseline 
with which measurements should be compared in order to 
see whether there are additional dynamical correlations. 
Clearly, this baseline can be improved by further impos- 
ing electric charge conservation. However, the number 
of charged particles in high energy heavy-ion collision 
is considerably larger than the number of baryons plus 
anti-baryons. Therefore, the corrections due to electric 
charge conservation will be sub-leading and only become 
relevant as the collision e nerg y is reduced. This will be 
discussed in detail in Ref. [26| . 

Ultimately it would be desirable to incorporate the ef- 
fects of global baryon number conservation into the vari- 
ous models or preferably into Lattice QCD. This is a very 



difficult task, and, therefore, we believe for the time being 
the present study will be helpful for the interpretation of 
the experimental data. 

In the next Section, we derive an analytical formula 
for the net baryon (and net proton) probability distri- 
bution under constraints imposed by global baryon num- 
ber conservation. Also the cumulant generating function, 
which can be used to compute cumulants of any order, 
is derived. In Section HI, we consider properties of the 
cumulants up to the sixth order and propose a new ob- 
servable, which is insensitive to the global conservation 
of baryon number. Comments and conclusions are pre- 
sented in Section IV and Section V, respectively. 



II. GLOBAL BARYON CONSERVATION 

Before we derive the relevant formulae let us remind 
ourselves what the problem at hand is. On the one hand 
we need to impose baryon number conservation on a sys- 
tem of baryons and anti-baryons following an underly- 
ing Poisson distribution. On the other hand we have to 
model the finite acceptance in an experiment, since for 
full acceptance the baryon number does not fluctuate. 
Here, we model the finite acceptances simply by a bi- 
nomial distribution, noting that in practice this may be 
more involved ^. 

These two tasks may be done in any order, i.e. one first 
derives the distribution for all particles subject to baryon 
number conservation and then folds with the binomial 
distribution or vice verse. Here we chose to first sepa- 
rate the system into observed and unobserved particles 
based on the binomial distribution. This is straightfor- 
ward since folding a Poisson distribution with a binomial 
results again in a Poisson distribution. Next we impose 
baryon number conservation on all particles, observed 
and unobserved. 

To get started, let us remind ourselves that the prob- 
ability distribution of the difference of two independent 
random variables, each drawn from a Poisson distribu- 
tion, is the so called Skellam distribution. Therefore, in 
our approach as well as in the HRG in the Boltzmann 
limit, in the absence of baryon number conservation the 
net baryon number is distributed according to the Skel- 
lam distribution (see e.g. [23, [23]) • Thus, in the following 
we will generalize the Skellam probability distribution by 
imposing global baryon number conservation. For the 
sake of simplicity, we perform our derivation for the net 
baryon number, and later we generalize our result to net 
protons. 

Suppose we have on average (Nb) baryons and {N§) 



^ Clearly, if all particles are observed, the baryon number will not 
fluctuate. 

^ For baryons and anti-baryons the classical (Boltzmann) approx- 
imation is justified due to a large ratio of mass to temperature. 



^ The use of a binomial distribution is correct if there are no cor- 
relations among baryons and anti-baryons other than the global 
conservation of baryon number which we take into account ex- 
plicitly. See section HVl for further discussion on this point. 



anti-baryons in the full phase space and an average net 
baryon number oi B = {Ng) — (Ng). According to our 
assumptions, both Nb and N^ follow a Poisson distri- 
bution. In order to model the finite acceptance we next 
split the full phase space into two subsystems, one repre- 
senting the measured particles, and one representing the 
unobserved rest of the particles. In the absence of any 
correlations particles are distributed between two sub- 
systems according to a binomial distribution, where the 
probability ps [Pb) ^o observe a baryon (antibaryon) is 
simply given by the fraction of the average number of 
observed baryons (antibaryons) to the average number 
of baryons (antibaryons) in the full phase space *. Con- 
sequently in both subsystems baryons and antibaryons 
are distributed according to Poisson distributions with 
appropriate means. 

As a result, the probability to observe ni net baryons 
in the measured phase space is again given by the Skellam 
distribution [27, 28] 
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and analogously in the unmeasured phase space, we have 



^2(^2) = M 



{1-Pb){Nb) 



"2/2 



{l-ps){Ns)_ 
xl„, (2zv/(1-ps)(1-Pb)) , (2) 

where A/1,2 are unimportant normalization constants and 

^ - V{Nb) {Ns). (3) 

The joint probability to have ni net baryons in the ob- 
served subsystem and 7^2 in the unobserved subsystem is 
given by 



P{ni,n2) = Pi{ni)P2{n2) 



(4) 



To impose the conservation of baryon number we multi- 
ply P{ni, 112) by 5ni+n'2.B and sum over all values of the 
unobserved net baryon number, 712, 

PB{ni)^NY Pl{ni)P2{n2)5n,+n,,B- (5) 

The normalization factor Af is fixed from the condition 
E PB{ni) = l. (6) 

Using Graf's addition formula ^ we obtain the net baryon 



probability distribution 
PB{n) -- 



^ /2 /, X (S-n)/2 

Pb\ f^-PB^ 



Pb 



(7) 



In {2z^pbPb) Ib-71 {^z^J{l -Pb){1-Pb)J 
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where z is given in Eq. (|3]). A detailed derivation of 
this result will be shown elsewhere [26|. The cumulant 
generating function g(f) for the cumulants Cfc 



f" 



^W^E^^IT 
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(8) 



is given by 



git) = ln(E„^B(")e"' 

'q+\ ^''^ lB{'2z^q+q_) 



= In 



(9) 



where q+ = I - pb + Pbs^ and q^ = I — pg +PBe^^- 

Finally, the probability distribution for net protons re- 
sults from the observation that all baryons other than 
protons (e.g. neutrons) may be considered as unobserved 
baryons. Thus, to obtain analogues of Eqs. ([7]) and (HJ 
for net protons, one simply defines the binomial proba- 
bility as follows 



PB = 



{nB) 
{Nb] 



{Ne 



(10) 



where (n^) is the mean number of observed protons (and 
analogously for pg) ^. 



III. CUMULANTS 

In this section we present the cumulants in the case 
where Pb = Pb — P- They result from Eq. ^ by taking 
the appropriate number of derivatives with respect to t 
and setting t to 0. 

As seen from Eq. ([9]), derivatives will generate Bessel 
functions of various order. Those can be simplified by us- 
ing the known properties of the Bessel functions ^ as well 
as the expression for the average number of baryons and 
anti-baryons in the case of baryon number conservation, 
{Nb)q and {Nb)(j, respectively, which will be discussed 
in the next section (see Eq. p^ *). For the first three even 



** Obviously < Pg jj < 1. 

^ Graf's addition formula 12911 is given by 
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which for t = 1 reduces to ^2 Ik(x)In—k{y) = In{x + y) 
k 



® Since there are at least as many neutrons in the system, 
(np>/(iVs> < i. 
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in B and do not depend on z. Thus, their ratios are 
uniquely defined in terms of p. 

Next let us define the ratio i?„ m as 



R„ 



(17) 



In Figs. 1 and 2, we show the dependence of the ratios 
of cumulants, Rn,m on p for realistic values of B and 
{■^B b)- The ratios of even cumulants are symmetric 
with respect to p —?' 1 — p as seen from Eqs. (|imi3l) . We 
note that for the Skellam distribution the ratios shown in 
Figs. 1 and 2 are unity. Therefore, we observe substantial 
corrections due to baryon number conservation. 

As already mentioned, the ratios of the odd order cu- 
mulants depend only on p. This allows us to construct 
the following combination 



D = R. 



i?3 



l--(l + 7)(3-7) 



(18) 



such that _D = for the baryon conservation corrected 
distribution Psin), Eq. ([T]), for any values of p, z and B. 

The upper (lower) sign should 
Also, D = for the 
Skellam distribution. Therefore, a deviation of D from 
zero may indicate physics, that is not related to global 
baryon conservation. 



Here, 7 = ±yr+ 8i?3,i- 

be taken for p < 3/4 {p > 3/4) ®. 



FIG. 1: Ratios of odd and even order cumulants as a function 
of the fraction of measured baryons, p. The parameters are 
B = 300, {Nb) = 400 and {Ng) = 100. 



order cumulants we obtain 

C2 = p{l^p){N)^, (11) 

C4 = C2 + 'S{p\^B^-cl) + 6pq{2z^pq~C2), (12) 

C6 = C4 +4(C4 - C2) - 10(2p(J + C2)(C4 - C2) 

-30pqip\^B^+cl), (13) 

where {N)^ = (A^s)c + (A^g)c, q = 1— pand z is defined 
in Eq. ^. Note that in the limit p -> one obtains 
cg « C4 w C2 ~ p{N)(j. It is worth mentioning that 
in case of the Skellam distribution, i.e., disregarding the 
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conservation law, we have c| 

The odd order cumulants can be expressed by polyno- 
mials in p 



(14) 
(15) 
(16) 

pB. As 
seen from Eqs. (1141116^ the odd order cumulants are linear 



ci = pB, 

C3 = ci(l-p)(l-2p), 

C5 = C3(l-12p(l-p)) 



For the Skellam distribution, cf 
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FIG. 2: Ratios of odd and even order cumulants as a function 
of the fraction of measured baryons, p in the range of values, 
which are of experiment interest. The parameters are B = 
300, (iVfl) = 400 and {Ng) = 100. 



* For an analysis of experimental data, the case with p < 3/4 
should be considered. 



IV. DISCUSSION AND COMMENTS 

Several comnients are in order regarding our results 
obtained in the previous sections: 



The distribution ([7]) depends on z = ^(A^b)(A^^), 
where (Nb) {{N§)) is the total baryon (antibaryon) 
number present in the Skellam distributions ([ij and 
©. Thus, (Nb) {{Nb)) is related to the system 
without baryon conservation. It is natural to ex- 
pect that baryon conservation will modify (Nb) 
{{Nb))-, however, as we argue below this correction 
is negligible. A straightforward calculation gives 



(Nb.b), 



I 



STl 



(2z) 



/b(2z) 



(19) 



where the upper (lower) sign corresponds to {Nb)c 
((iVB)c), with {Nb)c - {Nb)c = B. Here the 
subscript (•)p refers to averages obtained with full 
baryon number conservation. Under the constraint 
{Nb) — {Nb) = B, one can express z in terms of 
{Nb b)ci and to a very good approximation we find 



z^V{Nb)c-{Nb) 



BIC- 



(20) 



Using the properties of the modified Bessel func- 
tions, one can show that corrections to Eq. ((20)) 
are important only if both B and {Nb)c ■ {Nb)c 
simultaneously assume value of the order of one or 
smaller. This is never the case in heavy-ion colli- 
sions. Relation ([20l) together with the requirement 
that {Nb) - {Ng) = B ensures that {Nb) « {Nb)c 
and {Ng) ~ {Nb)c ^o very good precision. The 
same identities also hold if we only consider pro- 
tons. Therefore, the formalism developed in the 
previous section is of a great phenomenological 
value since it allows to calculate the effect of baryon 
number conservation on the probability distribu- 
tion and its cumulants given experimentally deter- 
mined average yields. This will be further elabo- 
rated in Ref. [26i. 

2. We have shown that the odd order cumulants do 
not depend on {Nbb)^ their ratios are indepen- 
dent of B and uniquely defined by one parameter, 
the fraction of observed baryons (protons) , p. This 
turns out to be very useful for the phenomenolog- 
ical analysis of experimental data. For example, 
chiral model calculations at non-zero baryon den- 
sities show that both R31 and R^^i are non triv- 
ial functions of temperature and chemical potential 
close to the crossover and the CEP. This is demon- 
strated in Fig. 3, where we present the results ob- 
tained in the Polyakov loop-extended Quark- Meson 
model [I4I for R^^i and i?5,i. We also show the new 
observable D, see Eq. ^TE\\ . which exhibits strong, 
temperature-dependent deviations from the base- 
line of £) = 0, even for temperatures below the 



pseudo-critical one, T < Tpc. Therefore, effects due 
to a possible phase transition should be accessible 
in experiment via an analysis of this new observ- 
able D. 
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FIG. 3: The ratios Ra,!, -Rs,! and D as a function of tempera- 
ture in the PQM model. The calculations are performed along 
the line of fixed hb/T ~ 0.5. Tpc is the crossover temperature 
at the corresponding chemical potential. 



3. For a given experiment, the parameter p can be 
roughly estimated. For example for the STAR ex- 
periment [llil the mean number of accepted protons 
in the most central Au-Au collision at ^/s = 200 
GeV in the measured phase space is approximately 
{up) « 7. The total mean number of protons can 
be estimated from data on dNp/dy at zero rapid- 
ity dNp/dy = 35 [30], and assuming flat rapidity 
distribution in the range —3 < y < 3. Therefore, 
{Np) « 35 • 6 = 210. We are, however, interested in 
the total baryon number, {Nb)- Therefore, {Np) 
is to be multiplied by some factor / that takes 
into account contribution of neutrons, A and other 
long living resonances. We estimated this factor in 
the thermal model. At the temperature T « 166 
MeV corresponding to y/s = 200 GeV we obtain 
/ = 2.5, so that {Nb) « 210 • 2.5 = 525. A similar 
number for {Nb) can be obtained using BRAHMS 
data [3l[ . The fraction of measured protons to the 
total number of baryons is p ss 7/525 ss 0.013. For 
this value of p and B = 350, we obtain i?4_2 « 0.95 
and i?6,2 ~ 0.77. There is some uncertainty related 
to the problem that we should only include those 
baryons that play a role in the quasi-equilibrium 
physics. This number is, however, difficult to esti- 
mate reliably. 

4. At sufficiently low energies {Nb) ^ {Nb) and 
z '^ B, the cumulant generating function ^ re- 



duces to ^ 



~g{t)^B\n[l-p{l-e*)] 



(21) 



In this case all ratios of the cumulants depend only 
on the fraction of observed baryons (protons) p. 
Taking (up) = 9 and (Nb) = 350 (number of par- 
ticipants in central collisions) we obtain p w 0.026 
and in consequence -R4.2 ~ 0.85 and i?6,2 ~ 0.32. 

5. In this paper we have disregarded other conserva- 
tion laws, e.g. electric charge conservation, which 
is expected to play a significant role at low energies. 
We believe, that at a collision energy higher than 10 
GeV, energy, momentum and electric charge con- 
servation can be neglected owing to high abundance 
of pions. A detailed investigation of electric charge 
conservation will be reported elsewhere 26] . 

6. We also have neglected non-equilibrium effects, ef- 
fects of interactions, volume fluctuations [32j re- 
sulting, e.g., from centrality fluctuations, etc. For 
example, here we have modeled the acceptance of 
baryons and anti-baryons simply by a binomial dis- 
tribution. This is only correct if there are no corre- 
lations among baryons and anti-baryon other than 
the global baryon number conservation, which we 
have accounted for explicitly in our calculations. In 
reality, one could very well imagine that the stop- 
ping of the baryons from the colliding nuclei may 
give rise to correlations and thus fluctuations of 
the baryon number which are not accounted for by 
the binomial distribution used in our calculation. 
These additional fluctuations will contribute to the 
various cumulants and thus need to be understood 
in order to extract any signal for critical fluctua- 
tions. Consequently, our proposed observable, D, 
will al so de viate from zero, since D is designed to 

9 For 2x < VkTT, IkCi^) ~ x''/r{k + 1). 



remove only correlations due to global baryon num- 
ber conservation. This in turn may be utilized to 
study baryon and anti-baryon correlations and fluc- 
tuations due to stopping. 

7. In the present paper we assume that baryon num- 
ber is conserved globally. However, it is plausi- 
ble that baryon number is conserved locally (a few 
units of rapidity). This effect would reduce the ef- 
fective total number of (anti)baryons (iV^ ^) and 
increase the value of pg ^ and, consequently, the 
corrections due to baryon number conservation. 



V. CONCLUSIONS 

We have studied the effects of global baryon conser- 
vation on the cumulants of net baryon and net proton 
fluctuations. We showed that the cumulants are substan- 
tially suppressed if global baryon conservation is taken 
into account. We also proposed a new observable that 
is insensitive to global baryon conservation but changes 
rapidly at the critical end point or the crossover. 
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